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Since  Peshkov^  has  confirmed  experimentally  Che  existence  of  a 

2 

temperature  wave  called  Che  second  sound  in  bulk  liquid  helium,  Landau 
evaluated  at  first  the  ordinary  second  sound  velocity  in  terms  of  the 
macroscopic  thermal  quantities  and  later  derived  the  second  sound  as  a 
collective  density  wave  in  the  elementary  excitation  spectrtm.  Transport 
coefficients,  i.e.,  thermal  conductivity,  viscosity  and  other  kinetic 

coefficients,  were  also  calculated  above  0.6  K  by  means  of  kinetic  equations 

3  4 

by  Landau  and  Khalatnikov. 

It  is  well  known  Chat  at  low  temperatures  and  low  pressures  the  thermal 
properties  of  superfluid  helium  are  dominated  by  low-momentum  phonons,  which 
do  not  possess  the  normal  dispersion,  but  rather  an  anomalous  dispersion.  In 
fact,  these  two  cases  are  quite  different  in  microscopic  processes.  The 
former  is  governed  by  four-phonon  processes  (4PP)  while  the  latter  is  by 

5 

three-phonon  processes  (3PP).  Recently,  several  authors  have  evaluated  the 
temperature  variations  of  Che  first  and  second  sound  by  using  the  wrong  normal 
dispersion.  However,  in  recent  articles,^  we  have  derived  a  Landau- type 
elementary  excitation  spectxrum,  which  has  anomalous  phonon-like  behavior  and 
rocon-llke  behavior  at  low  and  high  momenta,  respectively,  through  the  ring- 
diagram  approximation  in  two-  and  three-dimensional  liquid  helium.  In  the 
bulk  case,  the  obtained  excitation  spectrum  is  given  by 

E(p)  -  c^p[l  +  +  ...1 

E(P)  -  A  +  , 


3 


where  A,  and  ft  are  the  roton  parameters,  Is  the  sound  velocity  at 

absolute  zero  temperature,  and  5^  and  are  positive  constants  which  can  be 

determined  by  the  potential  parameters.  Here,  we  have  adopted  a  soft 

potential  with  a  Lennard-Jones- type  tail,  which  helps  make  a  smooth  connection 

between  the  attractive  part  and  a  soft  repulsive  core.  Using  these  excitation 

spectrums,  we  have  successfully  derived  and  explained  the  various  sounds,^ 

8  9 

sound  attenuations  and  transport  coefficients  in  two-dimensional  and  three- 

4 

dimensional  liquid  He. 

Concerning  the  coefficients  of  the  second  viscosity  in  bulk  liquid 

helium,  recently  Kirkpatrick  and  Dorfman^^  evaluated  these  cofficients  at  very 

low  temperature  (naA^  »  1)  on  the  basis  of  their  kinetic  equations^^  for  a 

dilute  superfluid,  where  n,  a  and  A  represent  the  number  density,  s-wave 

scattering  length  and  thermal  wave  length,  respectively.  In  a  previous 

12 

paper  wo  have  successfully  evaluated  the  coefficients  f2»  ^3  f4 
the  second  viscosity  in  two-dimensional  liquid  ^He  as  a  function  of 
temperature  by  solving  the  superfluid  hydrodynamics  through  the  theory  of  the 
kinetic  phenomena  developed  by  Khalatnikov. 

In  this  paper,  adopting  the  above  theory,  we  carry  out  a  three- 
dimensional  calculation  to  evaluate  the  coefficients  of  the  second  viscosity 
of  bulk  liquid  helium  via  3FF  in  the  region  of  very  low  temperatures ,  which 
was  not  done  by  Landau  and  Khalatnikov.  Near  zero  temperature,  the  important 
mechanism  is  3FP:  the  direct  process  of  emission  of  a  phonon  p  by  Pj^  ^2  ^ 
p  and  the  reverse  process  of  absorption  of  a  phonon  p  by  P2  +  p  ->  p^^.  The 
differential  decay  rate  for  3FF  in  three -dimensions  is  defined  as 


.  2s. 

“  w 


<F  Hjl> 


S(E^-E^) 


V  dpgdp 


(3) 
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where  is  Che  Hamiltonian  for  3PP.  The  transition  amplitude  between  the 
temperature -dependent  initial  state  1 1>  and  final  state  |f>  for  Che  direct 
process  is  given  by^ 


<F|H3|I> 


IL 

2  (2v)V2 


3c 


JL  pa 

2  dp  [p 


(n  . 

I  P  J  J  ?!  P2  P 


(4) 


Here  n^  is  the  phonon  distribution  function  with  momentum  p.  In  the  case  of 

the  reverse  process,  the  last  bracket  in  Eq.  (4)  should  be  replaced  by 
'  h 

[ (n  +l)n  n  ]  .  Then  the  total  decay  race  for  the  direct  process  becomes 

Pi  Po  P 


(5) 


where  u  -  (p/c)  (dc/d p)  is  Che  Grdneisen  constant.  For  the  reverse  process, 
n^  (n^  +l)(np+l)  is  replaced  by  (n^  ” 


Pi  P2 


Pi  P2  P 


From  Che  total  decay  rate  for  Che  direct  and  reverse  processes,  Che 
collision  integral  for  the  3PP  becomes 


irc 


3pp  2}lp 


2  r  _  _  *^^2 

(u+1)  PiPoP  -n  1  ' 

J  ^  ^  Plo  ^20  '  ^ 


(6) 


where  6n  is  equal  to  n  -n  ,  and  n  and  n  represent  the  equilibrium 


P  P 


Plo  P2o 
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distribucion  functions  for  the  phonons  with  monenta  and  p^.  respectively. 

Since  3PP  occur  largely  in  the  small  angle  events,  the  phonon  distribution, 
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which  depends  on  the  chemical  potential  a' ,  can  be  %n:itten  as 

n-  {exp[(o'+pc)/kgT]  -  1)’^  ,  (7) 

and  expanding  Sn  in  a  cower  series  in  terms  of  a'  ,  we  obtain  £n  -  n  •  n  -■ 

o 

4 

-n^(n^-fl)a’ .  The  relation  between  the  collision  integral  and  the  kinetic 
coefficient  can  be  expressed  by 


To  obtain  the  coefficients  of  the  second  viscosity,  we  follow  the 
calculation  in  the  third  way  given  by  Khalatnikov  and  show  below  only  the 
results  near  absolute  zero  temperature: 


(10) 


(11) 


6 


f3<T) 


'''■3PP 


U  ^£1: 

U  3fij 


(12) 


r^CT)  -  f^(T) 


(13) 


Here,  the  Grdneisen  constant  comes  from  differentiation  of  with  respect  to 

fi  and  can  be  taken  as  constant  near  absolute  zero  temperature.  rj^(T)  and 

r^(T)  are  equal  and  thus  satisfy  Onsager's  reciprocity  symmetric  principle. 

H^^(X)  is  the  number  of  phonons  per  unit  volume  at  temperature  T,  and  talejng 
4 

up  to  order  p  in  Eq.  (1),  we  can  express  this  as 


(14) 


where 


.3  Ihc  J 


“l 


3»6!r(6) 

4lf(4) 


“3 


5«8!r(8) 

4fr(4) 


(15) 


For  numerical  calculations,  we  have  adopted  the  parameters  which  vere 

deduced  from  the  analysis  of  the  elementary  excitation  of  the  bulk  liquid 
6  14 

helium  data,  ’  as  listed  in  Table  1.  We  have  also  taken  the  numerical  value 
of  u  to  be  1.8.^^  Using  these  parameters  and  substituting  Eqs.  (9)  and  (14) 


into  Eqs.  (10) -(12),  we  obtain 
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fj^(T)  -  X  IQ  ^ 

r2(T)  -  jj  ^  A(T) 

TjCT)  -  X  10  ^ 

where 


(16) 

(17) 

(18) 


A(T)  -  1  +  2a^T^  -  la^T^  +  (203+0^)!^ 
-  20^02!^  +  (a2+2a^a^)T:^  +  ... 


(19) 


in  Eqs.  (16)-(18),  all  fotir  coefficients  include  various  T-dependent  terms. 

The  T'^  dependence  also  appears  as  the  main  contribution  in  the  coefficient 

9 

of  the  first  viscosity, 


n(T)  - 


3!f(3)r(2)pc^  1 
Ag(u+1)2  T 


1.90  X  10*^  I 


(20) 


This  means  that  as  the  temperature  approaches  absolute  zero,  the  roton 
contribution  disappears,  and  thus  the  main  contribution  comes  from  3PP. 
Therefore,  the  temperature  dependence  of  the  first  viscosity  coefficient  is 
changed  dramatically  from  T'^  to  T'^. 

Kirkpatrick  and  Dorfman^®  have  introduced  the  classical  kinetic  theory 
technique  to  derive  two- fluid  hydrodynamics  in  an  expansion  of  the  unknown 
functions  in  terms  of  orthogonal  polynominals,  keeping  only  the  first  term, 
which  represents  the  first  Enskog  approximation  in  the  classical  theory. 

They  have  also  introduced  an  isotropic  linear  collision  operator  L  that  can 
be  divided  into  two  types  of  collision  processes;  an  excited  particle 
scatters  with  a  condensed  particle  to  produce  two  excited  particle,  or  vice 
versa  (^12^ ’  excited  particles  scatter  with  each  other  and  produce 

two  new  excited  particles  (122^ '  l^^se  two  processes  correspond  exactly  to 
the  processes  3PP  (p^,  P2  +  P3)  and  4PP  (p^^  +  P2  P3  +  P4)  •  The  collision 


8 


operator  L-^2  Isolds  only  for  very  small  momenta.  When  the  number  of  particles 
decreases  with  decreasing  temperature,  there  are  not  enough  collisions  to 
keep  a  local  equilibrium,  and  thus  their  hydrodynamic  solution  breaks  down, 
which  may  also  occur  in  very  dilute  classical  gases.  This  means  that 
Kirkpatrick  and  Dorfman  have  not  taken  into  account  3PP,  but  only  4PP. 

A  finite  lifetime  of  an  elementairy  excitation  in  a  quantum  liquid 
originates  from  two  processes:  one  is  the  collision  between  elementary 
excitations,  and  the  other  is  the  sponatneous  decay  of  a  phonon  into  two, 
three  or  more  phonons.  Since  the  collision  probability  tends  to  zero  at  very 
low  temperatures,  the  collision  process  becomes  less  important  and 
negligible.  At  very  low  temperatures,  the  main  excitations  in  liquid  ^He  are 
phonons.  On  the  other  hand,  we  note  that  3PP  are  the  most  important  process 
in  the  anomalous  elementary  excitation  spectrum  of  liquid  *^He  at  very  low 
temperatures  and  very  small  momenta. 

17  18 

Before  Kirkpatrick  and  Dorfman,  Ma  and  Popov  obtained  the  same 

2-9 

results,  i.e.,  a  damping  factor  on  the  order  of  q  T  ,  where  q  is  the  wave 

number  of  the  disturbance  in  the  calculation  of  the  hydrodynamic  eigenvalue 

describing  second  sound  propagation.  Comparing  this  with  that  of  Kirkpatrick 

and  Dorfman,  we  find  the  relation,  q  V/P^  -  q  T  ,  and  thus  have  rj  -  T  , 

where  is  proportional  to  T^.  Popov  has  also  derived  a  T  ^  dependence  of  rj 

through  the  functional  integral  analysis,  which  turns  out  to  be  mainly  the 

19 

same  contribution  from  4PP  in  Landau- Khalatnikov  theory. 

In  the  calculation  for  tj(T)  and  f(T)  given  by  Kirkpatrick  and  Dorfman, 
r;(T)  is  proportional  to  T‘^,  while  the  thermal  conductivity  k(T)  is 
proportional  to  T  ,  and  all  four  coefficients  of  the  second  viscosity  have 
only  a  T  dependence  for  very  low  temperature.  However,  Eqs.  (16) -(18) 
include  ''ot  only  T  terms  but  also  an  important  T  dependence.  As  T  0, 
their  four  coefficients  tend  directly  to  zero.  However,  Eqs.  (16) -(18)  show 
a  divergence,  which  is  quite  contrary  to  their  results.  We  emphasize  again 
the  difference  between  contributions  from  3PP  and  4PP.  The  discrepancy  in 
the  results  obtained  by  us  and  by  Kirkpatrick  and  Dorfman  is  due  to  the  fact 
that  their  hydrodynamic  solution  accommodates  only  4PP,  which  does  not  appear 
at  very  low  temperatures  in  the  Landau -Khalatnikov  theory,  and  breaks  down 
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when  3PP  are  Included,  while  our  theory  includes  only  3PP,  which  are  the  main 

contribution  to  the  transport  coefficients  near  zero  temperature  and  at  very 

small  momenta.  We  note  that  there  are  generally  thirteen  independent 

20 

dissipation  coefficients.  Here  we  only  consider  five  coefficients  and  do 

21 

not  discuss  the  pressure  and  frequency  variation  of  the  coefficients  of  the 
second  viscosity. 

Figure  1  illustrates  the  temperature  variations  and  magnitudes  of  the 
four  coefficients  of  the  second  viscosity  and  the  coefficient  of  the  first 
viscosity.  The  above  five  coefficients  are  all  positive.  For  very  low 

o 

temperatures  (well  below  T  <  0.6  K) ,  the  condition  <  C2f3  holds. 

Comparing  the  order  of  magnitudes  for  the  four  coefficients,  we  can  write 

^3  ^  ^1  ^  ^2’  four  coefficients  have  larger  values  than  that  of  the 

coefficient  of  the  first  viscosity.  We  note  that  the  thermal  conductivity 

(k)  and  five  coefficients,  as  mentioned  above,  for  thin  liquid  heliian  films 

12 

have  the  same  properties  as  those  for  bulk  liquid  helium. 

Figure  2  represents  the  values  of  the  coefficients'  of  the  first  and 
second  viscosity  as  a  function  of  temperature  in  the  functional  form  by  which 
they  appear  in  the  absoirption  coefficient  of  the  second  sound  given  by 
Khalatnikov.  The  behavior  of  this  curve  increases  moderately  as  temperature 
decreases  in  the  range  T  >  0.2  K  and  rapidly  diverges  near  zero  temperanire. 

It  shows  a  similar  behavior  to  that  in  the  high- temperature  region 
(T  >  0.6  K) . 

In  conclusion,  the  coefficients  of  the  second  viscosity  in  bulk  liquid 
helium  exhibit  a  T'^  dependence  for  temperaures  below  0.6  K,  like  that  of  the 
first  viscosity,  which  is  due  to  3PP  originating  from  the  anomalous 
excitation  spectrum  of  bulk  liquid  helium  at  low  momenta  and  low 
temperatures.  However,  there  are  no  experimental  results.  A  precise 
measurement  of  the  temperature  variation  of  the  second  viscosities  would  be 
desirable . 

This  work  was  supported  by  a  grant  to  Korea  University  from  the  BSRI, 
Ministry  of  Education,  1990,  Republic  of  Korea,  the  Center  for  Thermal  and 
Statistical  Physics,  and  the  US  Office  of  Naval  Research. 
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Table  1.  Parameters  for  bulk  liquid  helium 


n(A~^)  A/kg(K)  q^(A'^)  n _ c^(m/s)  _ gg(A^)  S^(k^) 

2.18x10'^  8.616  1.930  0.153m„  238.21  1.51±0.13  3.25+0.20  -16 
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Figure  Captions 

Figure  1.  Temperature  variation  of  the  coefficient  of  the  first  viscosity  (ij) 

and  the  four  coefficients  of  the  second  viscosity,  ^2 

fj*  The  curves  represent  only  the  magnitudes  of  these 

2 

coefficients,  where  the  units  are  g/cm»s  for  17  and  cm  /s  for 
(-f^)  and  cmVg'S  for 

Figure  2.  Coefficients  of  the  first  and  second  viscosity  as  a  function  of 
temperature  in  the  functional  form  by  which  they  appear  in  the 
absorption  coefficient  of  second  sound  given  by  Khalatnikov .  The 
quantities  and  are  the  superfluid  and  normal  densities, 
respectively. 


